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SOME LIMITATIONS ON ION-CYCLOTRON WAVE GENERATION AND 
SUBSEQUENT ION HEATING IN MAGNETIC BEACHES 
by Donald R. Sigman 
Lewis Research Center 

SUMMARY 

Thermal effects have been included in the analysis of ion-cyclotron waves generated 
by a Stix coil. Peak power coupling efficiency is reduced whenever ion or electron tem- 
peratures are sufficient for ion-cyclotron or Landau damping to occur. Because of the 
absorption of wave power under the coil, the wave power reaching an external "magnetic 
beach” is reduced. Other topics discussed include the effect of a Faraday shield on pow- 
er coupling and wave damping and interference effects which arise when waves are re- 
flected from the ends of the plasma column. 


INTRODUCTION 

In the ideal case, heating of plasma ions by ion-cyclotron waves can be viewed as a 
three -step process. First, a Stix coil (ref. 1, ch. 5) is used to couple radiofrequency 
energy to ion -cyclotron waves in a region of the plasma where the frequency of the waves 
(jo is less than the ion-cyclotron frequency Next, the waves propagate away from 
the Stix coil toward a region of lower magnetic field (called a magnetic beach), where 
the ion-cyclotron and wave frequencies are nearly equal. Finally, in the magnetic beach 
region, there are many ions which see the wave field at the cyclotron frequency and thus 
absorb energy from the wave. 

If the plasma ion and/or electron temperature is low underneath the Stix coil and in 
the immediate vicinity, collisionless damping processes (primarily electron Landau 
damping and ion-cyclotron damping) are not important; and the efficiency of coupling 
energy to the waves is very high (refs. 2 and 3) . However, if the plasma temperature 
is raised, there should be a point where electron Landau and/or ion -cyclotron damping 
will occur underneath the Stix coil. When this happens, it is certainly true that some of 
the ions will be heated outside of the magnetic beach region and not become trapped in a 
local mirror surrounding the beach (see fig. 1). It has also been Shown (ref. 4) that, 
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(a) Plasma - radiofrequency-coil geometry and finite- 
length system. 
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(b) Confining field with beaches. 

Figure 1. -Plasma model and direct-current confin- 
ing magnetic field. 


whenever there is damping underneath the coil, the coupling efficiency is changed. At 
Lewis this change in coupling efficiency has been seen experimentally in the Ion Cyclo- 
tron Resonance Apparatus (ICRA n). It was found that increasing the current in the Stix 
coil (and supposedly raising the temperature of the plasma) did not cause the power 
coupled to the plasma waves to increase at a rate proportional to the square of the cur- 
rent. Also, an energy balance study of the beach region of ICRA H gives experimental 
evidence that a significant amount of energy being coupled to the plasma does not reach 
the ions in the beach. 

This report examines the effects of finite plasma temperatures on coupling efficiency 
and determines the fraction of input energy reaching beach regions. Also, the effect of 
a Faraday shield on the wave damping is studied. Finally, a method is proposed where- 
by all the energy leaving the Stix coil can be made to flow in the same direction along the 
plasma column. 


SYMBOLS 


a c>ne-half Stix coil length 

B wave magnetic field 

Bq direct- current- confining magnetic field 
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wave electric field 
Stix coil current 

modified Bessel functions of the first kind 

Bessel functions of the first kind 

Stix coil current density 

modified Bessel functions of the second kind 

plasma wave number 

one -half system length 

electron density (equal to deuterium ion density) 

power coupled to plasma by Stix coil 

power flow from region of Stix coil 

plasma column radius 

energy loss rate from magnetic beaches 

radial position coordinate 

Stix coil radius 

electron temperature 

deuterium ion temperature 

volume of magnetic beaches 

thermal velocity parallel to confining field 

plasma dispersion function 
axial position coordinate 
axial center position of Stix coil 
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i argument of plasma dispersion function 

r energy loss time for beach regions 

£2 cu/J2. 

£2^ ion -cyclotron frequency 

co wave frequency 


METHOD OF ANALYSIS 

The plasma is assumed to be collisionless (ref. 2) and to consist only of deuterium 
ions and electrons. The particle density and the confining magnetic field are both as- 
sumed to be uniform. The ions and electrons are assumed to have Maxwellian velocity 
distributions, but the temperature of each species may be different. The plasma column 
is of finite length 2L. The waves are excited by a two -length Stix coil centered at z = z , 
where -L + (Xq/2) <z c <L - (Xq/2) and Xq is the wavelength of the Stix coil. When a 
magnetic beach exists, a negligible amount of the wave energy leaving the Stix coil is 
reflected back toward the coil and the plasma column appears infinitely long. This con- 
dition in the model is simulated by making L much greater than the characteristic 
damping length of the ion-cyclotron wave (L = 10 m is usually sufficient). Unless other- 
wise stated, L = 10 meters and z £ = 0. However, the effects of placing a reflecting grid 
near one or both ends of the Stix coil is also discussed briefly. 

The wave fields in this finite -length plasma column are represented as a Fourier 
series in a technique described in more detail in appendix A. There, the wave field 
equations are presented for two cases: a Faraday shield located at the plasma boundary, 
and no Faraday shield. 


RESULTS AND DISCUSSION 

Effect of Finite Ion and Electron Temperature on Coupling Efficiency 

Figures 2 and 3 give examples of the effect of finite ion and electron temperatures 
on wave coupling efficiency. The values of many of the parameters (which are listed on 
the figures) are typical of those in ICRA II. In figure 2 (where the electron temperature 
is held fixed at 10 eV), it is seen that coupling efficiency is reduced more than a factor 
of 3 by increasing the temperature from that of a cold plasma (T^ =10 eV) to that of a 
hot (but experimentally realizable) plasma (T^ = 1000 eV). 
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Figure Z - Power t ran sfer to plasma as function of ion- 
cyclotron frequency with ion temperature as a param- 
eter. Electron density, 2xl0 12 cm" 3 ; electron tem- 
perature, 10 eV; wave frequency, 40.8x10^ radians 
per second; coif current, 100 amperes; coil radius, 

7. 5 centimeters. 



Figure 3. - Power transfer to plasma as function of 
ion -cyclotron frequency with electron tempera- 
ture as a parameter. Electron density, 2xl0 12 
cm -3 ; deuterium ion temperature, 0.01eV;wave 
frequency, 40.8xl0 6 radians per second; coil 
current, 100 amperes; coil wavelength, 41 cent- 
imeters; coil radius, 7. 5 centimeters; plasma 
radius, 4 centimeters; deuterium gas. 
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It is expected that ion thermal effects become important when = (co - S2.)/kv ^ 

A ^ 

is of order one or less. Electron thermal effects are important when a q = co/kv^ is 

i ® 

of order one or less. The respective values of a_^ (at peak coupling for T. equal to 
10, 100, and 1000 eY are 3. 7, 1.8, and 0.8. For the electrons, respective values of 
Q!q for T g equal to 10 and 50 eV are 1. 41 and 0. 64. 

Figure 4 shows peak coupling efficiency as a function of density (again with T. as a 
parameter). At low densities the reduction in coupling due to increased ion temperature 



Figure 4. - Peak power transfer to plasma as function of electron 
density with ion temperature as a parameter. Electron tempera- 
ture, 10eV; wave frequency, 40.8xl0 6 radians per second; coil 
current, 100 amperes; coil wavelength, 41 centimeters; coil 
radius, 7.5 centimeters; plasma radius, 4 centimeters; deuter- 
ium gas. 


is large. This is because peak coupling occurs for small values of | o> - R. | . At higher 
densities s the peak coupling occurs at magnetic fields for which j cu - oj is larger and 
temperature effects are then less important. 

In general^ whenever the wave damping length is of the order of ? or less than ? the 
Stix coil lengthy a significant reduction in coupling efficiency occurs. 


Energy Absorption Underneath Stix Coil 

As was pointed out in the previous section, coupling efficiency is reduced when con- 
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ditions under the Stix coil are right for either electron Landau damping (co/kv^ ~ 1) or 

6 

ion-cyclotron damping ((co - fi.)AVn ~ 1). Obviously, whenever there is damping under 

the coil, not all the energy coupled -in leaves the coil region or gets to a magnetic beach 
located some distance from the coil. The amount of wave energy flowing away from the 
coil region has been calculated by evaluating the Poynting vector (axial component) at the 
ends of the Stix coil and integrating over the cross-sectional area of the plasma. The 
energy absorbed under the Stix coil is then the difference between the energy coupled-in 
P in and this integral of the Poynting vector P QU ^. The results of these calculations 
for some typical cases are shown in table I. 

Rows 1 and 2 of the table show that operating the Stix coil only slightly further from 
SI = 1 (0. 936 as opposed to 0. 952) can significantly increase the fraction of the input 
energy which leaves the coil (41. 5 percent as opposed to 17. 5 percent). 

When a Faraday shield is placed at the plasma surface, the wave E z field is forced 
to zero there: and subsequently the average E„ field across the whole plasma column is 
lower. When E z is less, the wave energy transferred to the electrons through Landau 
damping is less. This means that not only is the coupling efficiency increased, but also 
a higher percentage of the energy coupled -in leaves the region of the Stix coil in the form 
of waves. (Compare rows 2 and 3 and rows 8 and 9. ) 

Figure 5 shows the Poynting vector and wave E_ and E,, fields for the conditions of 

^ -j 

rows 2 and 3. Without the Faraday shield the energy flow ~(EXB)J is peaked near the 
radial boundary, with the peak moving slightly inward as z increases. Placing the 
Faraday shield at the plasma radius forces both E r and E z to be smaller near the 
surface and causes the energy flow to peak more towards the center of the plasma column. 

Row 5 shows that at fi = 1 both the coupling efficiency (P^ n = 3. 55 kW) and the frac- 
tion of the input energy leaving the coil (12. 3 percent) are low. The frequency of the 

7-1 

radiofrequency generator in ICRA II is presently 6. 5 megahertz (to = 4. 08X10 sec ). A 
comparison of rows 1 and 4 shows the improvement in performance that might be attained 
by increasing the frequency roughly 50 percent to 10 megahertz (w = 6.28x10 sec ). 

Row 7 gives results of a calculation at conditions similar to those for a fusion reac- 
tor. Increasing the density and increasing the temperature both have the effect of re- 
ducing the coupling; but this is offset by the necessity of operating further from Q = 1 
(because of high density) and by the use of a larger value of w (because of the larger 
confining field of the fusion plasma). 

All the calculations presented in figures 1 to 4 and in rows 1 to 7 of table I were made 
with the Stix coil at the center of a plasma column of length 2L = 20 meters. With 
L = 10 meters the waves for the conditions studied are damped out significantly by the 
time they propagate out to z = L; and thus there are no reflected waves returning to the 
Stix coil. When either or both of the ends of the system (which are assumed to be per- 
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Wave field E r , V/ cm Wave field E z , V/ cm Poynting vector, (EXB ) z , arb. units 



Radius, r/p 


Figure 5. - Poynting vector and wave E z and E r fields as 
function of radius. Electron density, 2xl0 12 cm" 3 ; 
deuterium ion temperatures, 100 eV; electron tempera- 
ture, 10 eV; ratio of wave frequency to ion -cyclotron 
frequency, 0.936; axial position coordinate, 42 centi- 
meters (unless otherwise noted). 



Ion-cyclotron frequency, rad/sec 

Figure 6. - Power coupling as function of ion- 
cyclotron frequency with Stix coil center 
position as a parameter. Electron density, 
2xl0 12 cm -3 ; deuterium ion temperature, 
100 eV; electron temperature, 10 eV; coil 
current, 100 amperes; coil radius, 7.5 
centimeters; plasma radius, 4 centimeters; 
wave frequency, 40.8x10° radians per sec- 
ond; Stix coil wavelength, 41 centimeters; 
one-half system length, 1000 centimeters; 
Faraday shield at r*p. 


9 



fectly reflecting) are brought close to the ends of the Stix coil, some of the wave energy 

is reflected back underneath the coil; and interference effects play a role in the coupling 

process. The result is that the coupling efficiency can be either increased or decreased 

depending on the position of the ends of the system relative to that of the coil. Figure 6 

gives power coupling as a function of SL, with L held at 10 meters and z (the position 

7 c _i 

of the center of the Stix coil) as a variable parameter. At = 4. 37X10 sec the coupl- 
ing is always higher than for the case of z = 0. The percentage of the energy leaving 

V 

the Stix coil does not change much (compare rows 8 and 2 or rows 9 and 3 of table I), but 
it is now all directed out one end of the coil. Thus it seems possible that such a reflector 



Figure 7. - Power coupling as function of ion -cyclotron frequency 
with system length as a parameter. Electron density, 2xI0 i2 
cm - 3j deuterium ion temperature, 100 eV; electron temperature, 
10 eV; coil current, 100 amperes; plasma radius, 4 centimeters; 
wave frequency, 40.8X10 6 radians per second; Stix coil wave- 
length, 41 centimeters; axial center position of Stix coil, 0; no 
Faraday shield. 
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at one end of the coil might cause more than twice as much energy to leave the other end 
of the coil, to be absorbed in a magnetic beach. 

Figure 7 gives power coupling as a function of £2^ with z Q = 0 and L as a variable 
parameter. It is seen that as L becomes less, the interference effects become more 
pronounced, and there are values of £2. where the coupling is much higher than in the 
case where L = 10 meters. It does not appear that this increased coupling could be taken 
advantage of in conjunction with a magnetic beach; but it does appear that it could be used 
for direct heating in the vicinity of the Stix coil. 


Energy Balance of Magnetic Beaches 

In the magnetic beach regions of ICRA II average ion temperatures of about 300 eV 

12 -3 

have been achieved with an electron density of 2x10 cm . The volume of both beaches 

3 -4 

is approximately 6000 cm . The measured energy loss time r is about 10 second. 

The energy loss rate, then, is approximately 

n (kT.)V 

R = — — « 6 kilowatts 

T 

This compares to an experimental energy input rate at the Stix coil of 26 kilowatts. 

Row 2 of table I is the best-known approximation to this experiment and predicts that 
41, 5 percent of the input energy leaves the Stix coil. (41. 5 percent of 26 kW is about 
11 kW. ) These two values (6 and 11 kW) are consistent with each other when it is con- 
sidered that (1) not all the energy leaving the coil gets to the beaches (which start about 
20 cm from the end of the Stix coil) and (2) the parallel ion temperature used for the 
theoretical calculation is not known to within better than a factor of 2. 

CONCLUSIONS 

A study has been made of the effect of finite ion and electron temperature on ion- 
cyclotron waves generated by a Stix coil. Whenever there is strong electron- Landau 
damping or ion-cyclotron damping underneath the coil, coupling efficiencies are re- 
duced. In a typical plasma (n = 2x10^ cm’^, T. = 100 eV, T = 10 eV, and 
co = 4. 08x10 sec ) only about 40 percent of the power coupled to the plasma column 
leaves the region of the Stix coil in the form of ion- cyclotron waves. Thus, the attainable 
ion temperatures in small local mirrors surrounding magnetic beaches is limited. For 
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plasma conditions close to those of a fusion plasma, coupling efficiencies can still be 
appreciable; although little power leaves the coil region. 

Reflections of waves from the ends of the system cause resonant cavity effects which 
produce greatly increased coupling efficiencies for certain values of fi.. If the Stix coil 
is near one end of the system, it may be possible to use reflections from that end to aid 
in increasing coupling efficiency and to direct a higher percentage of the wave power 
toward a beach region near the other end. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 14, 1971, 

129-02. 
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APPENDIX A 


FOURIER SERIES REPRESENTATION OF ST1X COIL 
FIELDS IN A FINITE-LENGTH PLASMA COLUMN 


Assume a uniform plasma column (length 2L, radius p) immersed in a uniform 
axial magnetic field Bq and surrounded by a two -wavelength Stix coil (length = 2 Aq, 
radius = s) centered at z = 0 (see fig. 1). Also, assume the plasma column is bounded 
on the ends by perfectly conducting plates so that the electric fields (E r , Eg) tangential 
to the plates may be set equal to zero. Consider the Eg component of the wave field 
generated by the coil (all other field components can be generated from Maxwell’s equa- 
tions once Eg is known). The current distribution in the coil j* is asymmetric and can 
thus be represented by a Fourier sine series. The field Eg is likewise asymmetric and 
can also be written as a sine series. Because of superposition, each term in the Eg 
series can be thought of as the field related to the corresponding term in the current dis- 
tribution series. The proper sine series to satisfy Eg = 0 at z = ±L is 


and 


Eg(r,z,t) = Eg^(r) sin^H^e ia,t 

n=0 


j*(r, z > t) 



j*6(r - s) sin 
n 



ia>t 

e 


Notice that, in order to satisfy the boundary condition at z = ±L for all values of r, it 

is necessary to make each term in the series equal zero. This is because the radial 

dependence of each of the coefficients Eg (r) is different. 

y n 

The Fourier coefficients Eg (r) have been derived previously as a function of j* 
(ref. 2). n n 
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where 


= [ a i J i(^i r ) + a 2 J i( y 2 r )_ 

k 

(r <P) 

= [ViM + ¥A r ) 

k 

(p < r < s) 

= c ^iW 


(r >s) 


_ -B ±Vb 2 - 4C 

%2= 2 

B _ k n (P + S) _ a) 2 p [ o> 2 (D 2 - S 2 ) 

2 c 2 c 2 S 

C = £k4 -2^k 2 S + ^(S 2 -D 2 ) 
s n c 2 n ^ 

and where S, -iD, and P are the hot- plasma dielectric tensor elements K^, K xy 
and K zz (ref. 1, ch. 9) (all functions of k n ). The computation of the dielectric tensor 
required evaluation of the plasma dispersion function Z(4). The approximations for 
Z(£) used in this report are found in appendix B. The coefficients a^, ag, bp bg, and 
Cg can be calculated by first computing the remaining E and B field components from 
Maxwell’s equations and applying the following boundary conditions: 

(1) Tangential components of both E and B are continuous across the plasma- 

vacuum interface (r = p) . 

(2) Tangential components of E and B are continuous at the external current sheet 

(r = s), except B z which is discontinuous by (47 t/c)]*. 

(3) When a Faraday shield is assumed to exist at the plasma boundary (r = p), E z is 

set equal to zero there and everywhere outside the plasma column. 

With no Faraday shield at r = p 5 
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a 1= 4l^5 Kl (ks) 


•'aPtyV* *1®®) + 5Kj(kp)] - kpljfKjP) I„(kp) - 5K 0 (kp)] 
y 2 P J i(^i p ) J 0 (^ 2 p ) " v iP 3 q(v 1 p)J 1 (v 2 p) 


a 0 = K 1 

^ 2 1 


K-i(ks) 


a = -4fficos K (kg) ^i^o(^lP) [ z i( k P) + «Ki( k P)] - kpJ^jP) [l Q (kp) - 6K Q (kp)] 

2 c 2 1 l/ 2P J l^lP^ J 0^2 p) " i; lP J 0^l p)J 1^2 p) 

^1 [r^PJo^jP) + pj 1 (^ 1 p)]j 1 (^ 2 p) - a 2 [y^PJp^^P) + PJ i(^ 2 P)] J l (y lP) 

k P K 0 (kp)K 1 (kp)r a 

r a = a l + PJ i(^P)] [y y 2 pJ 0^2 p ^ + J 1^2 P ^] 

- a 2 t /J/ 2 P J 0^2 p ^ + PJ 1^2 P ^] [y^PJ 0 ^l p) + J l^l p ^] 

K (kp) 

y = 

kpK Q (kp) 

S(k 2 + v *) - ^ (S 2 - D 2 ) 


S(k 2 + t/ 2 ) - fi. (S 2 - D 2 ) 


With a Faraday shield at r = p, 


_ 47Ti0JS K ’i( ks ) Q! 2 i; 2P J 0^2P^ 


“2 |/ 2 J o(‘'2 I>, 
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r b= t J l ( V )J 0 ( ^ a 2 y 2 P “ J 1^2^ J 0^1 p)a lV] k ^0 (kp ^ 


+ [j 0 (yiP)J 0 (^2 p ) I 'l pi '2 p ( Q! 2 " a i)] K i( k P) 

b 1 =[kpK 0 (kp)J 1 (i/ 1 p) + l / 1 pK 1 (kp)J 0 (y 1 p)Ja 1 + [kpK^kpJJ^p) + ^ 2 pK 1 {kp)J Q ( l / 2 p)ja 2 
b 2 = [ k P I 0 ( k P) J i( I, iP) - y 1 Pl 1 (kp)J 0 (^i p )] a x + [kpI 0 (kp)J 1 (i/ 2 p) - y 2 pI 1 (kp)J 0 (i^ 2 p)]a 2 


I t (ks) 

C n = " b -> + b n 

& -tr /l.„\ i Z 


K 1 (ks) 


The average power transferred to the plasma is given by 


-i/' 


E 0 (s, z) • j(z)(2?rs)dz 


Since j(z) is zero for L > z > a, we can write 


-if 

2 A 


E 0 (s, z) • j(z)(2ffs)dz 


And since E 0 (r, z) and j(z) are expanded in the same orthogonal series, it is possible 
to write 


P = -ITS 




Eg (s, z) • j n (z)dz 
n 


In references 2 and 3 the wave fields were calculated for an infinitely long plasma 
column, and thermal effects were neglected. In this case the fields were represented 
as a Fourier integral 


E(r, z, t) = 



E(r,k)e i < kz - < ‘ ,1 > dk 
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This integral was evaluated by a contour integral technique that gave solutions which 
were sums of the residues of the integrand. To evaluate the residues, the zeros of the 
functions T (see preceding page) had to be found; and then the derivatives of r with 
respect to k, dF/dk, evaluated. For cold plasmas, all the terms in r(k) are real and 
it is relatively easy to find numerically the values of k (all of which are real) which 
satisfy F(k) = 0. These k values are the wave numbers for the natural modes of the 
plasma column. When damping mechanisms (thermal effects) are included in the analy- 
sis, all the terms in F(k) may become complex; and in general the k values of the 
natural modes are complex numbers. It has been our experience that standard numeri- 
cal techniques do not work well in finding the roots of F(k) = 0 when damping is great 
and the real and imaginary parts are comparable. This is not to say it is impossible to 
find roots, just that it can be very time consuming. 

The advantage of the technique used in this report, then, is that the zeros of F(k) 
do not have to be found. For long plasma columns, the Fourier series solutions effec- 
tively amount to a numerical integration of the above Fourier integral. In addition, the 
effects of boundaries near the ends of the Stix coil can be included to the extent that one 
believes in the boundary conditions imposed. For power coupling calculations in which 
end boundary conditions are not important (such as when magnetic beaches exist), the 
technique is most applicable. 
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APPENDIX B 


APPROXIMATIONS OF PLASMA DISPERSION FUNCTION 


The calculations made in this report require the evaluation of the complex plasma 
dispersion function Z(4) with 4 being pure real. For smaller values of the argu- 
ment (4 < 3) and for the calculations presented in table I and figures 5 to 7, the continued 
fraction approximation of Derfler and Simonen (ref. 5) was used 


Z'(4) = -2iVff ? exp(-£ 1 2 ) - — £ 

1 + 1 




„ 1 e 2 
n - — £ 

2 


2n - — + 2n + — - 


Z(Q = lM±l 
24 

The curves in figures 2 to 4 are from a previous report in which, for small arguments, 
the two -pole approximation (ref. 6) was used. 


Z(S) = 




1 1 

q(o! - 4) o;*(a* + 4)J 

1. 10 


-= 0. 55 + 
a 



The change to the continued fraction was made because it was discovered to be more 
accurate in the range of interest. For larger arguments (4 > 3), and for all the calcula- 
tionSj the asymptotic expansion was used (ref. 1, ch. 8). 


Z(4) = iVn~ exp(-4 2 ) 


1 

4 
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